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We use a novel temperature-flow renormalization group technique to analyze magnetic and super- 
conducting instabilities in the two-dimensional t-t' Hubbard model for particle densities close to the 
van Hove filling as a function of the next-nearest neighbor hopping t' . In the one- loop flow at the van 
Hove filling, the characteristic temperature for the flow to strong coupling is suppressed drastically 
around t' c ~ — 0.33t, suggesting a quantum critical point between d-wave pairing at moderate t' > t' c 
and ferromagnetism for t' < t' c . Upon increasing the particle density in the latter regime the leading 
instability occurs in the triplet pairing channel. 
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In recent years the two-dimensional (2D) Hubbard 
model has mainly served as an idealized model in the 
efforts to understand the Mott-insulating, antiferromag- 
netic (AF), superconducting (SC) and more exotic prop- 
erties of the high-T c cuprates, like the pseudo-gap or 
stripe phases. While the latter phenomena are still under 
current debate and the Mott-insulating state is generi- 
cally restricted to larger onsite repulsion U, the AF and 
d-wave SC ground states in the Hubbard model can be 
understood, at least qualitatively, by calculations at weak 
or moderate U. Here, besides diagrammatic approaches 
taking into account subclasses of diagrams, like ladder 
summations [|l] and the fluctuation-exchange approxima- 
tion 0, one-loop renormalization group (RG) methods 
restricted to the saddle point regions ]||| and more re- 
cently taking into account the whole Fermi surface (^-0] 
have become a popular method to study the interplay be- 
tween magnetism and superconductivity in the Hubbard 
model. They indicate that for small nearest-neighbor 
hopping t' and close to half-filling, the weakly coupled 
state is unstable with respect to AF fluctuations, while 
at smaller or larger electron densities the RG flow ex- 
hibits a Cooper instability in the 4^2 _ y 2 -wave channel 
§. For t' -0.25t to -0.3t on the hole-doped side, 
the change from AF to d-wave instability leads over an 
intermediate regime, the saddle point regime of Ref. 

Originally the Hubbard model was introduced by Hub- 
bard, Gutzwiller, and Kanamori with a somewhat differ- 
ent motivation, namely as a model to describe the metal- 
lic ferromagnetism in the transition metals jl^,[ll| . Since 
then several attempts have been undertaken to clarify the 
conditions for itinerant ferromagnetism within the Hub- 
bard model and related systems. For several classes of 
one-dimensional |ll[] and higher-dimensional |l^] mod- 
els, spin polarized ground states can be proven to ex- 
ist, and also in the limit of infinite dimensions a number 
of clear results exist For the 2D t-t' Hubbard 

model at large U Becca and Sorella |l4|] recently found 



numerical evidence for a Nagaoka phase close to half- 
filling. At weak to moderate U however the situation is 
more involved. Hartree-Fock treatments fl5|| and more 
recently T-matrix approximation calculations combined 
with Quantum Monte Carlo |ll| and parquet approaches 
Jl7t indicated that for larger absolute values of the next- 
nearest neighbor hopping t' the ground state should be 
ferromagnetic (FM) in a certain density range around the 
van Hove-filling where the Fermi surface (FS) is near the 
saddle points of the band dispersion. A one-loop scaling 
analysis restricted to the saddle-point regions indicated 
similar tendencies [||. Nonetheless, one-loop RG calcula- 
tions |3]-f7| derived from exact RG equations and cover- 
ing the full FS have so far not indicated any FM ground 
states. Moreover the competition between FM and SC 
ground states remains an open problem. 

In this Letter we reinvestigate possible instabilities of 
the 2D t-t' Hubbard model towards magnetic and super- 
conducting phases from the weak-coupling perspective. 
The microscopic Hamiltonian we study is 

H = -* C M C J> ~ l ' c ls c i,» + U n ^ Ui i 

n.n., s n.n.n., s i 

with onsite repulsion U and hopping amplitudes t and t' 
between nearest neighbors (n.n.) and next-nearest neigh- 
bors (n.n.n.) on the 2D square lattice. We apply a novel 
one-loop RG scheme covering the full FS which we call 
the temperature-flow (T-flow) RG scheme. Like the ap- 
proaches in Refs. |^-|^] it is derived from an exact RG 
equation. However, no infrared cutoff is used, and in- 
stead the temperature T itself is the flow parameter. This 
new feature allows for the first time an unbiased com- 
parison between AF and FM tendencies, whereas RG 
schemes with a flowing infrared cutoff artificially sup- 
press particle-hole excitations with low wavevectors. A 
detailed discussion of this issue and a derivation of the 
T-flow RG equations are given in Ref. |J. Using the 
RG formalism, we obtain a hierarchy of differential equa- 
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tions for the one-particle irreducible n-point vertex func- 

(n) 

tions as functions of the temperature. Integration 
of this system gives the T-flow. We pose the initial 
condition that at a higher temperature To the single- 
particle Green's function of the system is described by 
Go(iuj,k) = [iu> — e(fc)] -1 and the interaction vertex is 
given by a local repulsion U. This is justified if T is 
large enough. We truncate the infinite system of equa- 
tions by dropping all with n > A. In the present 
treatment we also neglect selfenergy corrections and the 
frequency dependence of the vertex functions. This re- 
stricts the scheme to a one-loop equation for the spin- 
rotation invariant four-point vertex 1^ . Starting with 

weak to moderate interactions, we follow the flow of 1^ 
as T decreases. 

T^ is determined by a coupling function Vr(fci,fc2, 
fea) For the numerical treatment, we define elon- 

gated phase space patches around lines leading from the 
origin of the BZ to the (±7r, ±7r)-points, and approxi- 
mate Vr{ki, &2, /C3) by a constant for all wave vectors 
in the same patch. We calculate the RG flow for the 
subset of interaction vertices with one wavevector on 
the FS representative for each patch, with initial con- 
dition Vr (ki, k,2, £3) = U. Most calculations were per- 
formed using 48 patches. Further we calculate several 
static susceptibilities, as described in Ref. (7). Here we 
focus on the most divergent ones, namely d 2 .2„ y 2-wave 
and p-wave pairing susceptibilities, the AF susceptibility 
Xs{q = (tt, tt)) and the FM susceptibility Xs(Q —> 0). This 
way we analyze which classes of coupling functions and 
which susceptibilities become important at low T . For a 
large parameter range, we observe a flow to strong cou- 
pling. The approximations mentioned above fail when 
the couplings get too large. Therefore we stop the flow 
when the largest coupling exceeds a high value larger 
than the bandwidth, e.g. Vr, max = 18£. This defines a 
characteristic temperature T c of the flow to strong cou- 
pling. Because breaking of continuous symmetries is im- 
possible in 2D at T > 0, we interpret T c as an estimate 
for the temperature where ordering occurs when a small 
additional coupling in the third spatial direction is in- 
cluded. 

Here we describe the results for the RG flow of the 
coupling function for the 2D Hubbard model with initial 
interaction U — it and varying value for the next-nearest 
neighbor hopping t' . The chemical potential is fixed Jll| 
at the van Hove (VH) value fi = At' so that the FS always 
contains the saddle points at (±7r, 0) and (0, ±7r). 

We decrease t' from to more and more negative val- 
ues. The overall behavior of the flow to strong coupling 
is shown in Fig. |l|. Focusing on AF, d-wave and FM sus- 
ceptibilities, we observe three distinct parameter regions 
with typical FS shown in Fig. ||. The first regime is closer 
to half-filling, where t' > — 0.2t. Here T c is relatively high 
and the AF susceptibility grows most strongly (see Fig. 



. The dominant scattering processes are AF processes 
between FS parts connected by wavevectors w (tt,tt). 
When we further decrease t', T c drops continuously, and 
for t' < —0.2t, the cZ-wave susceptibility takes over as the 
leading susceptibility. For t' rss —0.25t and band fillings 
slightly larger than the VH filling we again find a regime 
where the flow of d-wave SC and AF processes is strongly 
coupled. This is the saddle point regime studied in Ref. 

. Here we focus on the flow when we go to even more 
negative t and at the same time adjust the chemical po- 
tential such that the FS stays at the saddle points. Then 
T c drops by several orders of magnitude for t' < — Q.3t, 
while for t' < — 0.33i it rises again. For these values 
the flow to strong coupling is dominated by processes 
with small momentum transfer and the FM susceptibil- 
ity Xs(Q = 0) ls by far the most divergent susceptibility 
at low temperatures (right plot in Fig. ||; note that the 
FM susceptibility dominates already at a higher T, where 
the largest couplings are still well within the bandwidth. 
The overall behavior strongly suggests a quantum critical 
point around t' c « — 0.33t where the ground state changes 
from d-wave singlet SC to FM. 
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FIG. 1. The characteristic temperature for the flow to 
strong coupling, T c , versus next nearest neighbor hopping am- 
plitude t' from a 48-patch calculation (the crosses show data 
for 96 patches). The chemical potential is fixed at the van 
Hove value /1 = At' . T c is defined as temperature where the 
couplings reach values larger than 18i. As criterion for the dis- 
tinction between antiferromagnetic, d ;c 2„ y 2-wave pairing and 
ferromagnetic regime we use the temperature derivatives of 
the susceptibilities at the scale where the couplings become 
larger than lOt. 

Lin and Hirsch Jll| analyzed the Hartree-Fock phase 
diagram for the t-t' Hubbard model and found the FM 
state to be stable against the AF state for t' < — 0.324i. 
Since they used the bare Hubbard repulsion, the d x 2_ y 2- 
wave phase is absent. The critical i'-value from Hartree- 
Fock is practically equivalent to the critical t' c « — 0.33t 
from the RG treatment. On the other hand in Hartree- 
Fock one would find a first order transition at t' c , while in 
the RG the characteristic temperature gets suppressed to 
lowest values, hinting at a quantum critical point. A T- 
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matrix approximation (TMA) applied by Hlubina et al. 

indicated a stable fully polarized FM state beyond 
a critical *' « —0.43*. The screening in the particle- 
particle channel of the TMA is included in our T-flow 
scheme, which gives a somewhat larger window for a FM 
state with nonzero polarization at the VH filling. The 
density range around the VH filling where FM dominates 
gets more and more narrow as *' is decreased towards t' c . 
Recently Irkhin et al. [jl7| applied a parquet scheme with 
a simplified dispersion and found indications for FM at 
the VH filling for very similar values of *'. 

A crucial feature of models where large-spin ground 
states have been found flo|-|l2|| is a sharp peak of the 
density of states (DOS) at the bottom of the band. This 
property is shared by the 2D t-t' Hubbard model for t' m 
—0.5*, but at the van Hove filling, the DOS is also peaked 
at the Fermi level. In our RG calculation the DOS at the 
Fermi level, and not that at the bottom of the band, 
is essential for the FM tendencies: they also exist, and 
dominate, for smaller |*'| in the range described above, 
and they are strongly reduced when the Fermi level is 
raised above the VH energy (see below) . This sensitivity 
of the FM regime to the FS location points to another 
effect which is not included in the present calculation. It 
is probable that the inclusion of selfenergy effects in the 
flow will shift the FS and alter the low energy DOS such 
that the actual FM regime might be shifted, or changed 
more strongly. The analysis of these effects is left for 
future work. 
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FIG. 2. Flow of the AF (dashed line) , d-wave pairing (solid 
line) and FM (dotted line) susceptibilities in the AF nesting 
regime (left plot, [i = — O.Olt), d-wave regime (middle plot, 
li = —t) and FM regime (right plot, fi = —l.St). The crosses 
(triangles) on the horizontal axes denote the temperatures 
where the largest couplings reach 5t (lOi). The insets show 
the Fermi surfaces and the 48 points in the Brillouin zone for 
which the coupling function is calculated. 

Next we investigate the flow to strong coupling for 
densities away from the VH filling. The RG flow for 
smaller |*'| is very similar to the flow found with the 
momentum- shell techniques which has been extensively 
discussed in Refs. Therefore we now focus on the 

case *' < — 0.33* corresponding to the FM regime at the 
VH filling and search for SC tendencies when the band 
filling is varied. In the 2D Hubbard model the occurrence 
of triplet pairing in vicinity of a FM phase has been pro- 



posed by several authors |]T^,|o). Our approach allows a 
more detailed analysis of this interplay. 

The Cooper pair scattering Vr{k,—k,k') transforms 
according to one of the 5 irreducible representations 
of the point group of the 2D square lattice. Among 
those there is the one-dimensional representation with 
d x 2_„2 -symmetry important for smaller \t'\ and the two- 
dimensional representation transforming like p x and p y . 
In a solution of the BCS gap equation, components be- 
longing to different representations do not mix, while e.g. 
within the p-representation components p x oc cos 9 and 
p y oc sin 9 and corresponding higher harmonics oc cos 39 
and oc sin 39 (also called /-wave) will occur together. 

As shown in the right plot of Fig. || for t' — — 0.45*, 
the critical temperature drops by several orders of magni- 
tude when we increase the particle density per site from 
the VH value (n) « 0.47 at \i = -1.8* to (n) w 0.58 
at (J, = —1.7t. Further upon moving away from the VH 
filling, the growth of the FM susceptibility gets cut off 
and the p-wave triplet SC susceptibilities with symme- 
try p x oc cos9 or p v oc sin# diverge at low temperature. 
Higher order harmonics oc cos 39 and sin 39 diverge in a 
weaker fashion. The pair scattering at low T in this pa- 
rameter range is shown in Fig.[|. One nicely observes that 
the pair scattering involving particles close to the saddle 
points is suppressed as the odd-parity nature of the p- 
wave pairing requires an opposite sign of the pair scat- 
tering Vr(k, —k, k') e.g. for k' (0,7r) and k! rs (0, — n). 
Comparing the temperature derivatives of FM and SC 
susceptibilities when the leading couplings become larger 
than 12*, the transition from the FM to the p-wave SC 
regime occurs at a particle density (n) ~ 0.52 per site or 
fx = —1.76*. Again, just as in the interplay between AF 
and d x 2_ y 2 SC fluctuations for small |*'|, our one-loop 
flow to strong coupling exhibits a smooth evolution from 
the FM dominated to the p-wave dominated instability. 
This suggests a transition between the two types of or- 
dered states as a function of the band filling provided 
that symmetry-breaking at T > becomes possible in 
a three-dimensional environment. From our analysis the 
transition appears to be first order, but additional in- 
teraction effects beyond our one-loop calculation could 
change that. 

With our method we cannot calculate the SC gap func- 
tion directly, but we can use the pair scattering obtained 
from the RG flow as a pair potential in the BCS gap equa- 
tion. In our case, it is plausible to assume |El| that the 
SC gap function will be a nodeless superposition of the 
two components with symmetries p x and p y , e.g. given by 
A(fc) = Ao(fe) (k x + ik y ), which maximizes the condensa- 
tion energy. The real-valued prefactor A (k) takes care 
of the anisotropy within symmetry-related FS parts. The 
direction of the Cooper pair spin remains indeterminate 
in the absence of spin-orbit coupling. The strongly angle- 
dependent gap function calculated with the rescaled pair 
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scattering taken from the RG flow is shown in Fig. ||. 
A similar SC order parameter is under debate for the 
superconducting state of Sr 2 Ru04 p2j. In this quasi- 
2D system one of the three Fermi surfaces has a similar 
shape as the case studied above. It will be interesting to 
investigate the RG flow and the SC properties suggested 
by that in a more realistic three-band model. 



dominated by triplet SC correlations suggesting a highly 
anisotropic energy gap. 
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FIG. 3. Left: Flow of ferromagnetic (solid lines) and 
p^-wave (dashed lines) susceptibilities for t' = — 0.45t at 
fi — —1.78t (thick lines) and [i = —1.744 (thin lines). Right: 
Characteristic temperature T c where the largest couplings 
reach 18t versus chemical potential fi. 
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FIG. 4. Left plot: Pair scattering V(k, —k, fe) with k vary- 
ing around the Fermi surface for t' — —0.454 and /i = —1.71 
at T = 1.1 ■ 10~ 6 t. For the solid line ks is fixed close to 
the saddle point (— 7r,0) while for the line with the dots, 
kg ~ (—1.06, —0.86) is near the zone diagonal. Middle plot: 
Real and imaginary part of the SC gap function obtained from 
the solution of the BCS gap equation with the rescaled pair 
scattering from the RG flow. Right plot: gap magnitude. 

In conclusion, using the new temperature-flow RG 
scheme, we have given a comprehensive analysis of the 
flow to strong coupling for the 2D t-t' repulsive Hubbard 
model on the square lattice. The flow to strong coupling 
close to half-filling for small next nearest neighbor hop- 
ping t' with its AF and d x 2_ y 2 -pairing regimes is qualita- 
tively similar to the RG flow known from the momentum- 
shell techniques . For larger absolute values of t' be- 
yond a quantum critical point at t' c w — 0.33t the flow to 
strong coupling is dominated by FM fluctuations. This 
has not been found with the momentum-shell schemes 
used previously for the reasons explained in Ref. 

||. Finally, for larger \t'\, when the electron density is 
increased slightly above the VH density, the FM tenden- 
cies get cut off at low T and the flow to strong coupling is 
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